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The three-nucleon scattering problem in a nuclear medium is considered within the Faddeev 
technique. In particular the deuteron break-up cross section that governs the formation and the 
break-up reactions of deuterons {NNN <-> Nd) in a nuclear environment is calculated at finite 
temperatures and densities. A significant enhancement of the in-medium break-up cross section 
with increasing density has been found. 
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Formation of light clusters such as deuterons, helium and alpha particles is an important aspect of heavy ion 
collisions at intermediate energies, see e.g. Empirical evidence, including recent experimental data on cluster 
formation, see refs. [^|J^, indicate that a large fraction of deuterons can be formed in heavy-ion collision of energies 
^ ■ E/A< 200 MeV. 

I During the expansion of the system the density can drop below the Mott-density of deuteron dissociation 

■ In this region the deuteron abundances will be determined by deuteron formation, NNN dN, and break-up, 
dN — > NNN, reactions. Since the deuteron formation rate can be expressed through the break-up probability, the 

, final outcome of the reaction will be essentially controlled by the deuteron break-up cross section. Previous studies 
^ of the kinetics of deuteron production have utilized the impulse approximation to calculate the reaction cross section 
OO , at energies above 200 MeV/A 0. For lower energies, viz. E/A < 200 MeV, the impulse approximation fails and 
CO • a full three-body treatment of the scattering problem is necessary. Furthermore, a consistent treatment of cluster 

I formation in expanding hot and dense matter requires the inclusion of medium effects into the respective reaction 

• cross sections. 

\^ [ The essentials of the three-body problem in the vacuum case are well known, see e.g. ref. In the following 
. we utilize the AGS formalism |^ suitably modified, in order to treat the three-body problem taking into account 
the nuclear medium in the quasi-particle approximation. To this end, we will rely on a separable representation 
of the nucleon-nucleon potential. This choice simplifies the problem considerably. A systematic investigation of 
separable parameterizations of "realistic" potentials has been pursued e.g. by Plessas and collaborators |^0|. We note 
O [ that solutions of the three-body problem using "realistic" NN potentials have been achieved e.g. by the Bochum 
2 ■ group ijll I , and the Bonn group in the framework of the M^-matrix approach . 

Recently, separable approximations have been used to solve the pion deuteron scattering problem and success- 
fully applied to coherent photo- and electro-production of pions on the deuteron [Q. A detailed discussion of the 
numerical method used in these cases has been given in [ p^ . In the present work we have modified this pion deuteron 
code to solve the nucleon-deuteron scattering problem. Our results for the dominant three-body transition matrix 

■ elements agree within ~ 5% with the previous calculations by Dolcschall . 

The AGS equations for the transition matrix Uaf3{z) are given in a compact notation by 

3 

Uii^{z) = (1 - 6c,f3)G^\z) + ^(1 - ^„^)T^(z)Go(z)C/^„(z), (1) 

7 = 1 

where the Greek indices represent particle channels. In the absence of the medium the free three-body resolvent is 
Gq{z) = {z — Hq)^^ with z = E + ie, where Hq is the Hamiltonian of the noninteracting three-body system and E 
is the three-body energy. The three-body transition matrix T^{z) in the channel 7 is related to the t-matrix of the 
two-body subspace Tj{z) (denoted by "hat") via, 

(q^p^|r^(2;)|qT,PT,) = {q^\q^){p^\f^{z^)\pj) , (2) 



*on leave of absence from Physics Department, Yerevan State University, Armenia 



1 



with = z — 3g^/(4m), where m is the nucleon mass, and (q-y) are the three-body momenta in the center of 
mass frame for the pair (odd) nucleon, respectively. For a separable two-body potential, which is given by (overall 
particle channel index 7 suppressed) 

Vr ^ \gr)K{9r\, (3) 

where r denotes a set of quantum numbers providing a complete specification of the two particle state, the two-body 
i-matrix is then given by 

Tr{z) = \gr)Tr{z){gr\. (4) 

Solving the two particle Lippmann-Schwinger equation determines -?> (see e.g. 

Tr{z) ^ (^X;' ~ {gr\Go{z)\gr)y\ (5) 
Multiplying the AGS equations with {gf3m\Go{z) [GQ{z)\ga,n)] from left [right] they reduce to [p^ 

l3n,am 



where 



{gf3n\Go{z)U0aGo{z)\gam), (7) 
Zf3n,am = (5/9ii|(l " Spa)Go{z)\gam) ■ (8) 

We restrict the two-body channels to the dominant ones, i.e. ^5*0 and "^Si — "^Di. For the separable ansatz we use 
the parameterization of Phillips jl^ . The parameters are taken from Ref . , which shows that gross features of the 
elastic and break-up cross sections as well as the differential elastic cross section up to Eiat = 50 MeV are sufhciently 
well reproduced. To calculate the break-up cross section we use the optical theorem. Figure |l| shows our results for 
the integrated cross section along with the experimental data . 

Note that the agreement between the present calculations on free nucleon-deuteron scattering and the experimental 
data can still be improved, through using more realistic potentials and including contributions of higher partial waves. 
However, our present goal is to study the general features related to the inclusion of nuclear medium, which are 
expected to be quite significant, cf. p2| . 

The influence of the nuclear medium on the three particle problem is treated in mean-field approximation. Here, 
two effects appear, i) self energy corrections and ii) Pauli blocking. The single particle propagator is obtained within 
the Hartree-Fock approximation of the self energy. The Hartree-Fock energy shift is given by {V{12, 12) = (12|y|12)) 

A^^(l) = [^(12, 12) - 1/(12, 21)] /(2), (9) 

2 

where the multi-indices 1 , 2 denote momentum and other quantum numbers of the respective particles, the Fermi 
function is given by /(ei) — {exp[(ei — fi)/T] + 1}~^, and the single particle energy by ei — kf/2m + A^^{1). The 
Hartree-Fock shift is evaluated using the same potential as for solving the AGS equations. To simplify the three-body 
treatment we introduce the effective mass approximation. The quasi-particle dispersion relation is then written as 
£i — kf/2m* + Aq^ , where the constant term A^-'^ can be absorbed into the chemical potential fi. The effective 
mass m* is density and temperature dependent, and decreases from the value m* /m = 1 in the zero density limit to 
e.g. m*/m = 0.97 at T = 10 MeV for the Mott density of the deuteron n = 8 x 10"^ fm"^. 

The Bethe-Salpeter equation for the three-body Green's function (7(123, 123; z) in mean field approximaton, in- 
cluding Pauli blocking and self energy shifts, is obtain from the Matsubara Green's function technique as 

G(123, 123; z) = ^ ~ ■^(^) ~ -^(2) ~ /(^) (123|123) 

z — ei — £2 — £3 

+ V G(123, 1'2'3'; z) ( V{l'2\ 12) (3'|3} ^ ~ -^^^^ ~ -^(^^ + cycl. perm.^ (10) 

lf^3, V Z- 61^62- £3 J 

This leads to modifications of the quantities entering in the AGS equations. Pauli blocking then appears in the 
transition matrix T^{z) and the resolvent Go{z). The thermodynamic f-matrix acting in channel 7, see eq. (0), has 
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been derived within the Matsubara-Green's function techniques, see Ref. It turns out that, in the low density 
approximation used here, the effect of the Pauli blocking may be absorbed into an effective potential. Inclusion of all 
medium corrections effectively results in the following replacements in the eqs. (||), (||), (^), and (H), viz. 

m —>■ m* , (11) 

dam ^ 9am ~ (Qa) ^ 9am- (12) 

Here (Q) is the angular averaged Pauli blocking factor, 

(Q(p, q)) - 2^ / ^ dx[l - /(eq/2+p) - /(eq/2-p)]. (13) 

where x = p • q/ (pq). The center of mass of the three nucleon system is set equal to the momentum of the medium. 
The AGS equations are properly symmetrized for identical particles and then solved in the rest system of the medium. 

We present our results with amplitudes normalized to the residue of the two-body ^-matrix in the deuteron channel 
d, i.e. 

N;j EE Res(f,z = E;^) = {9UGliE;f}\9;dr\ (14) 

where E*^ is the in-medium deuteron binding energy. The advantage is that the optical theorem, which is used to 
calculate the break-up cross section, has the same formal structure as in the vacuum case and is valid without any 
further modifications. For an alternative formulation of the in-medium optical theorem in the two-body case see e.g. 
Ref. jsj. The in-medium break up cross section a'^j,{^Eiab) that depends on density n and temperature T is then 
defined using the amplitudes N^^X^^ adi^)^ad resulting from solving eq. (^) with the replacements given in eqs. ( [Tl| ) 
and ([l^ ) in the same way as in the vacuum case, see e.g. 1^ or fl^ , p^ . 

Our results are depicted in Figure ^ The solid lines represent the free space break-up cross section as shown in 
Figure]^. The dashed lines show the break-up cross section, cr* rp(Eiab) for densities n = 0.1, 1, 3, 5, 7 x 10"'^ fm~^, 

respectively, as a function of the laboratory energy Eiab- The Mott transition occurs at the density n = 8 x 10"'^ 
fm-3. 

The medium effects significantly modify the vacuum break-up cross section. Two qualitative features are observed. 
First, the break-up threshold is shifted towards lower scattering energies with increasing density of the nuclear matter. 
This kinematical effect is due to the decrease of the deuteron binding energy with increasing density. For an illustration 
of the dependence of the threshold energy in the laboratory system i?*'* on different temperatures and densities see 
Figure |^. Second, the cross section increases considerably with increasing density. The maximum is enhanced by one 
order of magnitude for the largest density value considered. For densities larger than the Mott densities the deuteron 
disappears as a bound state. 

Also, it is instructive to see how the medium dependent cross section converges to the "free" cross section. At 100 
MeV the deviation of the in-medium cross section from the free in this model is in the order of 10%. From inspection 
of Figure || we conclude that the dominant changes in the cross section takes place at rather moderate energies, i.e. 
where the impulse approximation fails, and the Faddeev technique has to be used. 

To summarize, we have calculated the deuteron break-up cross section, cr^ j^{Eiab), using the AGS formalism 
modified in order to take into account the effects of quasi-particle energies and the phase space occupation effects. It 
has been shown that there are strong modifications of the deuteron break-up cross section compared to the vacuum 
scattering case. A study of the reactions dynamics within a kinetic approach to the cluster formation p5|-|27[|, which 
will include the deuteron break-up and formation probabilities based on the present work, would provide an improved 
description of the physics of deuteron production in heavy-ion reactions. 
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G.R. and A.S. are thankful to V.B. Belyaev and W. Plessas for discussions. 
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FIG. 1. A comparison of the total, elastic, and break-up cross sections with the experimental data of Ref. [21]. 



5 



1500 



1000 



Sit 







I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 

^ a: p=0.0001 fm'^ 

b: p=0.001 fm'^ : 
c: p=0.003 frn ^ 
d: p=0.005 frn ^ 
e. p =0.007 fm'^ I 

r \ : 

■ \ 

■ \ 
\ 

\ ^ 



500 - 








20 



40 



60 



80 



100 



EJMeV] 



FIG. 2. Break-up cross section at temperature T = 10 MeV. Free cross section is shown as solid line. Other lines are due 
to different nuclear densities, see text. 
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FIG. 3. Shift of threshold energy E in the laboratory system due to nuclear density n at different temperatures T. 
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